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Abstract 

We define a new family of matrix product states which are exact ground states 
of spin 1/2 Hamiltonians on one dimensional lattices. This class of Hamiltonians 
contain both Heisenberg and Dzyaloshinskii-Moriya interactions but at speci- 
fied and not arbitrary couplings. We also compute in closed forms the one and 
two-point functions and the explicit form of the ground state. The degeneracy 
structure of the ground state is also discussed. 
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1 Introduction 



In the past few years, a lot of interest has been attracted to the subject of matrix prod- 
uct states [HG2E1H1H1EJE1E1E1E]. ^ nere nas been a revival of interest in this subject 
[TUl [ITJ [12j [131 (13 05] due to the advances in quantum information theory [17] and 
the techniques developed in this field. Quite recently the method of matrix product 
states has been used for obtaining new models of many body systems for which an 
exact ground state can be defined and its many-body properties like the entanglement 
between various sites can be calculated in closed form [181 1151 1201 12H 1221 1231 1231 1251 12B]. 
This is interesting in view of the fact that in some works, transitions in entanglement 
behaviour has turned out [301 EQ to be a good signal of quantum phase transitions [29J . 

The reason for this high level of interest is the complementary role that the fields 
of condensed matter physics and quantum information play in investigation of many 
body systems. The matrix product formalism [32} [33] is one of the subjects which 
lies at the borderline of these two subjects. As is well known, in this formalism, 
one starts from proposed states whose expansion coefficients are the trace of prod- 
uct of given matrices. While for numerical investigations, i.e. the density matrix 
renormalization group (DMRG), one usually starts from large dimensional matri- 
ces, to simulate ground states of given Hamiltonians, in the approach which is used 
for finding exactly solvable models, one starts from low dimensional matrices and 
finds family of Hamiltonians for which these states are exact ground states. This 
is the approach which has been used in many of the works in the past few years 

In a recent work [26J, we classified all spin 1/2 matrix product states which arise 
from two-dimensional matrices and showed that some of them can undergo phase 
transitions. Doing the same kind of classification for higher dimensional matrices 
turns out to be a formidable task. A general classification in effect requires classifi- 
cation of solutions of a system of non-linear equations, which obviously is extremely 
difficult. However it is possible to define new models by starting with special classes 
of matrices which have nice algebraic properties. Such properties allow us to bypass 
the above-mentioned difficult problem. The price that we will pay is that we will 
obtain only restricted family of models and not all models. However if we start from 
a nice algebra, it is usually possible to end at physically interesting models and by 
the subsequent analysis of these models beyond the matrix product formalism, many 
useful information can be obtained for such models which are otherwise difficult to 
obtain. We will later see explicit examples of this method. 

Following this method, in this paper we will introduce one such family of matrix 
product states which is based on a simple algebra 

XZ = ujZX (1) 

where lv is a root of unity. We use a finite dimensional representation of this algebra 
in which Z and X are represented by generalized Pauli operators. Everything in 
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the model, from the state and its parent Hamiltonian (i.e. the Hamiltonian which 
has the matrix product state as its ground state) to the one- and two-point func- 
tions depend on the dimension of the representation. We will show that while for 
D = 2 dimensional representation of the algebra, the parent Hamiltonian, is a one 
parameter of Heisenberg spin chain, for D ^ 2 dimensional representation the Hamil- 
tonian contains both Heisenberg and Dzyaloshinskii-Moriya interactions [2T|, [28] , but 
for specified values of couplings which depend on the dimensional representation D. 

It is to be noted that the matrix product state is by construction a translation- 
invariant state and in some cases, it happens that this state is a sum of two ground 
states each of which breaks the translation symmetry. This is the case with Majumdar- 
Ghosh |34[ [35] model [20] . Here we will see a similar phenomena, where the matrix 
product state is the sum of many ground states each of which breaks the translation 
symmetry of the original Hamiltonian. 

The calculation of the one- and two-point functions is an important step in this 
analysis. While such a calculation is straightforward in the formalism of matrix prod- 
uct formalism, a close scrutiny of their structure will give clues as to the whole or 
part of the ground space structure and by following these clues we can obtain valuable 
information about the degeneracy structure of the ground state. 

The structure of this paper is as follows. In section ([2]) we give a brief introduc- 
tion to matrix product states, in section ([3]) we introduce the representation of the 
algebra ([1]) and construct the family of matrix product states and derive the parent 
Hamiltonian, in section (|55[) we calculate the one and two point functions. Finally 
we end up with a discussion. 

2 A brief introduction to matrix product states 

First let us review the basics of matrix product states. Consider a homogeneous ring 
of N sites, where each site describes a d— level state. The Hilbert space of each site 
is spanned by the basis vectors i = 0, • • • d — 1. A state 



is called a matrix product state if there exists D dimensional complex matrices A{ 6 




(2) 



C DxD , i = 0---d-l such that 




(3) 



where Z is a normalization constant given by 




(4) 
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and 

d-l 

E--J2 K ® U (5) 

i=Q 

Here we are restricting ourselves to translationally invariant states, by taking the 
matrices to be site- independent. 

Let O be any local operator acting on a single site. Then we can obtain the 
one-point function on site k of the chain {^/\0(k)\^) as follows: 

mo(km = «fA ,0) 

where 

d-l 

E := Y,(i\0\j)A*®A r (7) 

i,j=0 

The n-point functions can be obtained in a similar way. For example, the two-point 
function (^>\0(k)0(l)\^f) can be obtained as 

www = ^ggpgg) (8 ) 

where E {k) := E k - x E E~ k . Note that this is a formal notation which allows us 
to write the n-point functions in a uniform way, it does not require that E is an 
invertible matrix. Also by considering the permutation operator P defined as 

P | a(3) =| /3a) (9) 

and using equation ©we see that E* = PEP and since P 2 = P we find that Z* = Z, 
meaning that the normalization is real as it should be. Also for a Hermitian operator 
O in equation © the one point function (\I/|0(A;)|\E r ) will be real. The same reasoning 
applies for n-point functions. 

2.1 The Parent Hamiltonian 

Given a matrix product state, the reduced density matrix of k consecutive sites is 
given by 

tr((A* - --A* ® A h ■ ■ ■ An, )E N - k ) 
Piv~ik,h~h - tr(E N ) ' 

The null space of this reduced density matrix includes the solutions of the following 
system of equations 

d-l 

E c jl ... jk A jl ---A jk = 0. (11) 

Given that the matrices Ai are of size DxD, there are D 2 equations with d k unknowns. 
Since there can be at most D 2 independent equations, there are at least d k — D 2 
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solutions for this system of equations. Thus for the density matrix of k sites to have 
a null space it is sufficient that the following inequality holds 

d k > D 2 . (12) 

Let the null space of the reduced density matrix be spanned by the orthogonal vectors 
\e a ), (a = 1, • • • s, > d k — D 2 ). Then we can construct the local hamiltonian acting 
on k consecutive sites as 

s 

h:=Y^ Va\e a ){e a \, (13) 

a=l 

where fj, a 's are positive constants. These parameters together with the parameters of 
the vectors |e») inherited from those of the original matrices A%, determine the total 
number of coupling constants of the Hamiltonian. If we call the embedding of this 
local Hamiltonian into the sites I to I + k by hii+k then the full Hamiltonian on the 
chain is written as 

N 

H = J2Ki+k- (14) 
1=1 

The state \^f) is then a ground state of this hamiltonian with vanishing energy. The 
reason is as follows: 

N 

{*\H\V) = tr(H\y)(V\) = tr(h lt i +kPlt i +k ) = 0, (15) 

1=1 

where pi k+l is the reduced density matrix of sites I to l+k and in the last line we have 
used the fact that h is constructed from the null eigenvectors of p for k consecutive 
sites. Given that H is a positive operator, this proves the assertion. 

In view of the above introduction, we have a clear receipe for constructing matrix 
product states and a family of parent Hamiltonians. First one chooses the range 
of interaction k and then choose appropriate matrices, throwing away all spurious 
degrees of freedom by appropriate transformations A\ — > SAiS -1 and reducing 
further the degrees of freedom by imposing symmetries. In this way one ends with 
a reasonable set of matrix product states, which hopefully may have applications in 
description of real physical systems. In the following section, we want to study a 
particular family of such states arising from the simple algebra HJ 

3 The model 

In this paper we specify the auxiliary matrices Aq and A± as D dimensional gener- 
alization of spin (Pauli) operators in x and z directions. This is a D dimensional 
representation of the algebra [TJ Hereafter, for convenience, we denote Aq and A\ 
respectively by X and Z. In the so called computational basis for qudits, spanned by 
the orthonormal vectors {|0), |1), • • • \n)} and we have 

d-l 

X = Y,\n + l)(n\, (16) 

n=0 
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and 

D-l 



n=0 



2ni 

where cj = e o 



Z=J2" n \n)(n\, (17) 
One finds from (|16l and [T7|) that 

ZX = ujXZ. (18) 

As explained in section^ equation ([5]), the transfer matrix for our model has the 
following simple form 

E = X ®X + Z® Z. (19) 

The matrix E plays a central role for calculating the normalization of the ground 
state and also in determining all the correlation functions. We need to find the 
spectrum of this matrix. The right and left eigenvectors of E turn out to be 

1 D - 1 

\Mr)}-= - 7 =5>- rn Kn + fc>, (20) 
<V> fc (r)|:=^£V>,n + H (21) 

V D n=0 

both of which correspond to the eigenvalue 

A fc (r) =u k + u r , k,r = 0, 1, - D - 1. (22) 
It is easily checked that {ipk( r )\' t Pk , ( r ')) = ^k,k'^r,r'- 

Also we need to know Z = tr(E N ), the normalization constant which from (|22p 
turns out to be 

Z D (N)= (u r +u k ) N , (23) 

k,r=0 

where we have re-labeled the partition function as Zd{N) to emphasize its dependence 
on the number of sites N and the dimension of the representation D. From the 
binomial expansion, we find 

Zd{N) = E E f 1 ) ^ k(N - l) (24) 

r,k=0l=0 \ ) 

Summing over r and k and noting that Y,?=o U}Vl = D&l,mD and Ylk=o w k ( N ~ 1 ' = 
D5n-i m iD where m and m! are integers we find that Zjj(N) is nonzero only when N 
is a multiple of D. Hereafter we assume that this is the case, that is we set iV — > ND 
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in all the equations. The physical necessity of this requirement will become clear when 
we discuss the explicit form of the ground states. With this modification, we find the 

N n/v 

Z D (DN) = D 2 J2i Dl )• (25) 
1=0 

This is a sum of ./V terms, where N is the number of sites and as such is not a 
close expression. However we can use an identity and rewrite it as a sum of D terms 
which is much simpler. To this end we consider the following binomial identity 

ND ND 

(l + o/)^=£( N D )u (26) 

1=0 

and sum over both sides for s = 0, 1, • • • D — 1 to obtain 



D-l D-1ND 

I 

I 



s=0 1=0 

ND T. -r 7 — . D-l 



=0 s=0 1=0 



= U N F)Y," Sl ( 2 7) 

1=0 s=0 

However, since co D = 1, the sum over s in the right hand side is non- vanishing only 
when I is a multiple of D and in that case it will be equal to D. Therefore we find 
that 

D-l 

Z ND = Dj2(±+u S ) ND (28) 

i=0 

In contrast to (|25p . this is a sum over D terms which is independent of the system 
size and obviously is much simpler. 

For example for D = 2 and for D = 3 we have respectively 

Z 2N = 2 2N+l (29) 

and 

Z 3N = 3(2 m + 2(-l) N ) (30) 

We will now consider the parent Hamiltonian, for which the matrix product state 
is an exact ground state. 

3.1 The Parent Hamiltonian 

To find the parent Hamiltonian, we use the relation CijA{Aj = which in view of 
([!]) takes the form 
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C 00 X 2 + C iXZ + C W ZX + C n Z 2 = 0. (31) 

For O / 2 the operators X 2 and Z 2 are independent. However in any dimension 
ZX = uXZ and hence the vanishing of the left hand side gives the following relations 
on the coefficients 

Coo = C n = 0, Coi + w(7 10 = 0. (32) 

According to ([32]) , these equations define the null space of the density matrix of two 
adjacent sites. This null space is spanned by only one vector, which we denote by |e), 

|e) = |10) - w|01). (33) 

Thus the local Hamiltonian is given by 



h k ,k+i = tz(I ~ ^o-fc+i) - ua+a k+1 - ua k (35) 



h k ,k+l = l e X e U,JH-l = (I 10 ) " M01})((10| - 57(01|) fc)fc+ i. (34) 
Writing the above expression in terms of the Pauli operators, we find that 

1 

2 l 

After a simple re-scaling H — > 2H and neglecting an additive constant, the full 
Hamiltonian becomes 

N 2vr 2vr 
H = Y, cos -p( a M+i + " a k^k+i + sin — - <r y k a% +1 ). (36) 

k=l 

This is the general form of the Hamiltonian H when we use the D ^ 2 dimen- 
sional representation of the algebra ([T]). Therefore this Hamiltonian contains both 
Heisenberg and Dzyaloshinskii-Moriya interaction \27\ 128] . but with specific and not 
arbitrary couplings. 

In case we use the D = 2 dimensional representation, a different Hamiltonian will 
be obtained, since in this case we have X 2 = Z 2 = I and the null-space condition 
([3"7]) will have a different solution. In this null-space condition (]37]) turns out to be 

(Coo + C n )/ + (Coi - C W )ZX = 0, (37) 

and the null space is spanned by two vectors, namely 

|ei> = |00> — |11>, |e 2 ) = |01) + 1 10). (38) 

This larger null space means that the Hamiltonian has now two arbitrary coupling 
constants, and by a re-scaling H — > 2H, the Hamiltonian takes the form 

N 

H = J2(Jl - J2) Ha z k+1 - a%at +l ) + (J, + J 2 )(J + a\a\^). (39) 
i=i 

In this case the Dzyaloshinskii-Moriya interaction is absent. Note the Hamiltonian [39] 
cannot be obtained from [3H] in a certain limit. Up to now the parent Hamiltonian has 
been obtained for both types of representations of the matrix product algebra (j37p . 
The next step will be to calculate the correlation functions. We limit ourselves to one 
and two-point functions. 
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4 One and two-point correlation functions 



The matrix product formalism allows a straightforward calculation of correlation func- 
tions. According to ([7]), the matrix product operators corresponding to the spin 
observables a x , a y and a z are 



Ex 

Ey 

E K 



X <g> Z + Z ® X, 
-iX <g) Z + iZ (g> X, 
X ®X - Z ® Z. 



(40) 



For finding the average of Pauli operators in each site of the chain we use equation 
(jSJ). As explained in section (|2J), for finite systems, this requires a diagonalization of 
the transfer matrix E, while in the thermodynamic limit even this problem simplifies 
to the determination of the largest eigenvalue of the transfer matrix. In some cases, 
i.e. for the two-dimensional representation of the algebra (pQ), one can even proceed in 
an alternative and very simple way. Therefore we organize the content of this section 
in two separate subsections, where we deal separately, using different methods, with 
2-dimensional and higher dimensional representations of the matrix product algebra. 



4.1 Two-dimensional representation 

Let us start by noting that for the two-dimensional representation, we have X 2 = 
Z 2 = I and ZX = —XZ. The same relations hold with X — > Y. These relations and 
the expression for the transfer operator E given in (|19p leads to the simple relations 

E X E = E y E = E Z E = 0, (41) 

from which we immediately obtain that 

{af> = (4) = K> = 0- (42) 

It is obvious from the relations (|4ip that any correlation function containing a x 
and (j y operators on non- adjacent sites will also vanish. For adjacent sites we have 
to calculate traces of operators like E 2 E N ~ 2 , E 2 E N ~ 2 and E x E y E N ~ 2 . It is simply 
obtained from (j4"0j) that 



which leads to 

(of o-fc+i) = 



ElE = 0, 



E 2 E 



trE 



-4E, 



2N-2 



ElE = 0, 



trE 2N 



1, 



(43) 



(44) 



where we have used equation (|40j) . The last relation implies that the ground state 
is one in which the the spins of adjacent sites are aligned in the y directions but in 
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opposite senses. There are two such states and the matrix product state, which by 
construction should be translation invariant, is the sum of these two states, namely 

\*) = ^(\+y,- y r N + \-y,+yr N ), (45) 

where \± y ) := ( ^ are the two eigenstates of the a y operator. By using the 
relations 

a z \± y ) = \T y ), (T x \± y ) =±t|=F»>, v y \±y) = ±|±tf>, 

one can readily check that the two states in the superposition of are actually 
annihilated by the Hamiltonian, so they are the ground states of the Hamiltonian. In 
fact a simple way to see this, is to rotate all the spins in the odd-numbered sites by an 
angle ir around the z-axis. Under such a rotation, the Hamiltonian ([39]) transforms 
to 

N 

H' = £(Ji - J 2 ) {a z k al +l + a%at +l ) + {J x + J 2 )(/ - a^ +1 ). (46) 
i=i 

or equivalently by adding and subtracting terms of the form o v k o v k+v it transforms to 

TV 

H' = J2(Ji - J2) fa • ffk+i) - UivM+v (47) 
i=i 

where ■ (fk+i = a k a k+i + a k a k+i + a k a k+i anc ^ we have dropped a total additive 
constant. In case that < J% < J2, the first term is a ferromagnetic interaction which 
tends to align all the spins in one single direction, the direction is determined by the 
second y — y interaction. Hence we have a doubly degenerate ground state of the form 
I $+) := \y+)® 2N and |#_) := \y-)® 2N for H' . Rotating back all the odd-numbered 
spins by tt around the z axis gives the two ground states of H as expressed in (|45p . 
where their sum has been obtained as a matrix product state. Thus in this case the 
MPS formalism guides us through the complete degenerate structure of the ground 
state. However if J\ > J 2 > 0, then the first term in H' is an anti-ferromagnetic 
interaction with a highly degenerate complex structure. In this case the MPS only 
gives a very small part of the spectrum. 



4.2 Higher dimensional representations 

We now turn to higher dimensional representation where the Hamiltonian contains 
a Dzyialoshinski-Moriya interaction. In this case the commutation relations do not 
allow a simple calculation of the correlation functions for finite systems. Therefore 
we calculate the correlations directly in the thermodynamic limit. In this limit, only 
the largest eigenvalue and eigenvector of the transition operator E will survive the 
limit N — > 00 and one finds that 
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(Ok) = T (^max\Eo\^max), (48) 

■^max 

and 

(OfcOfc+l) = T2 (^max\Eo\\max) , (49) 

^max 

where X ma x and \X ma x) are the largest eignevalue of E and its corresponding 
eigenvector. From (|22p we know that the largest eigenvalue of E is 2 corresponding 
to the generalized Bell state |^oo) = 77^ Z)n=o For the one-point functions we 

need the matrix elements of the operators E x ,E y , and E z on this Bell state. It is 
obvious from (j4*0|) that (V'ool-Krl^oo) = (V'ool^lV'oo) = 0. For E z we find that 



D-l 

D 
1 



(ipoo\E z \i/j o) = — (^mlXigiX-ZfgiZln.n) 



n,m=0 
D-l 



Therefore we find 



jy ^2 (fim,n+l ~ $m,n) = 0. (50) 
n,m=0 



(a*) = (<?y) = {*z) = 0. (51) 



We now proceed to calculate the two-point functions for adjacent sites. To do 
this we need the following matrix elements which are readily obtained from (j40j) after 
simple manipulations: 



(^00 \E%\ tpoo) = uj + lj' 1 = 2 cos— , 
(ipoo\Ey\ipoo) = u + UJ' 1 = 2 cos— , 

(iP 00 \E 2 z \^ 00 ) = 0, ' (52) 



and 



, 2-7T 

(i/joo\E x Ey\i/joo) = i(u)-0J ) = -2sin— , 

1 2-7T 

(ij)oo\EyE x \il)Qo) = -i(uj-u) ) = 2sin— . (53) 

Moreover all the other quadratic matrix elements turn out to vanish. From the 
above matrix elements, one can easily obtain the correlation functions from (I49p as 
follows: 

1 2-7T 

KX+i> = = 2 cos ~jj > (^"M = 0, (54) 

and 



2-7T 

(^ + i) = -(^M + i) = 2 sin c- (55) 
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.A, .A. 
* / \ / \ 

Figure 1: (Color online) One of the product states \x) (for D=6). The matrix product 
state is a superposition of D such states, each shifted by one site with respect to the 
previous one. 



All the other two-point correlation functions vanish for nearest neighbor sites. The 
above results can be written in the following compact form: 

(a k ■ o-fc+i) = cos — {a k x a k+1 ) = -sm— z. (56) 

From these correlation functions along with the result (|45p which we obtained for 
the explicit form of the ground state for the two-dimensional representations, we can 
guess an explicit form for the ground state of the Hamiltonian ()36|) . We will later 
check the validity of our guess by direct calculations. The ground state (135]) . turned 
out to be a superposition of two product states, each one of them was nothing but a 
juxtaposition of spins in opposite directions. In fact each of the two product states are 
ground states of the Hamiltonian separately, and the matrix product state superposes 
them to ensure its own translational invariance. The same thing may happen here, 
that is we may find product states each of which is a ground state of the Hamiltonian 
and the matrix product state is a superposition of them. Inspired by the form of the 
ground states for the D=2 representation, we construct a product state which gives 
the same form of correlation functions as in (|56p . This construction is a generalization 
of the states in (135]) . Let 

\<j)) := s'm(j)\+ x ) + cos<j)\+y), 

where \ + x ) and \+y) are the spin states in the positive x and y directions respectively. 
We now form the following product state on D consecutive sites (instead of two sites 
in the D = 2 case) 

\xm-.= \m-^-^)--M- 2 -^^). (57) 

Each consecutive spin has been rotated anti-clockwise by an angle of around the 
z axis, in analogy with the D = 2 dimensional case. That is, the spins in each block 
of D consecutive sites are distributed regularly around a circle in the x — y plane. It 
is obvious that for any such block state, comprising D consecutive sites, the average 
of spin operator in any direction should vanish, that is 

1 D_1 

== „ E (x(^)WiW)) =0, a = x,y, z. (58) 
i=0 

The reason is nothing but the regular distribution of the spin states \<f>) over the 
circle. The same type of reasoning, which can be supported by direct and explicit 
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calculations, shows the validity of the two-point correlations (|56p for any block state 
of the above form. 



We can extend such a state to the whole lattice as a non-translation invariant 
Ix^))®^, which obviously has the same kinds of correlations as in (|56|) . The matrix 
product state, being translation invariant is nothing but the following superposition 

\*d) ■= \x(o)r N + \x{ 2 -^T N + \x{% )T N ■■■ + \x( 2{L } D 1)T - ))«* (59) 

This formula for the MPS ground state generalized the explicit formula (I45D to higher 
dimensional representations of the algebra. 

It is instructive to know the origin of the character of the ground state (|59p . Let 
us proceed as in the D = 2 case and rotate the spins individually, but this time, 
according to a different pattern. Consider the local Hamiltonian 



2tT 27T 

h k ,k+i = cos — {a%at +1 + a y k a y k+1 ) - a z k a z k+1 + sin — {a%a y k+l - <J y k at +1 ) (60) 

and rotate the spin at site k + 1 by an angle 9 = — jy around the z axis. Under this 
transformation, the spin operators transform as follows: 



a k+i — ► a k+l 
erf +1 — > - cos 
a y k+1 — > sin 6 o% +1 -cos Q (61) 



o-fe+l — > - cos 6 a% +l + sin 9 a y k+1 



Under this transformation, it is readily seen that h^^k+i transforms to an isotropic 
Heisenberg ferromagnet, h k ^fi 

ftgk = " + ^K+i + ° Z k°t + i) ■ (62) 

It is well known that the ground state of the Heisenberg Hamiltonian H Heis = 
J2k hkk+i 1S ^ ne st a t e in which all the spins align in one direction, and since the 
Hamiltonian has rotational symmetry, the ground state consists of a whole spin mul- 
tiplet, which is the multiplet with the Highest spin. For a lattice with ND sites, this 
multiplet has the spin whose top state is ^p-) = \z+)® ND . The other states 
are obtained by the action of the angular momentum operator L_ := Ylk=l a k on 
this top state. That is, the other un-normalized ground states are 

|™, m)=i ?^-»'|^,^). (63) 

Acting on these states by the inverse rotations fc i2 2 (^ L ) we obtain the ground 
states of the original Hamiltonian in the form 
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Ill 




Mf-m) ND ND 



2 ' 2 



(64) 



The ground states in |x(^))® from which the matrix product state (|59p is con- 
structed are only part of the ground space of the Hamiltonian. As expected in the 
matrix product formalism one starts from one single state, but a careful analysis of 
that single state and its correlation functions, assisted by some guesswork and physical 
analogies, can unravel the whole degeneracy structure of the Hamiltonian. 

5 Discussion 

We have started from a simple matrix product algebra, namely X Z = u Z X where 
uj d = 1 and proceeded to obtain the properties of the matrix product state which 
corresponds to this algebra. There is a distinct difference between the cases D = 2 
and D ^ 2. In the former case, the model describes a system of spins on a line 
with Heisenberg type interaction but with prescribed couplings, in the latter case, 
the model describes both a Heisenberg and a Dzyaloshinski-Morya (DM) interaction. 
In both cases we have been able to calculate the one and two point functions and from 
the insight that these functions have provided, we have been able to obtain a more 
thorough understanding of these models, including the degeneracy structure of their 
ground state. One possible line of expansion of this work is to start from a Heisen- 
berg ferromagnetic chain of arbitrary spins in a magnetic field. If now we rotate in a 
periodic pattern the magnetic fields and demand, similar to the work of Kurman et al 
[36J, that all the cross terms except the Heisenberg and the DM interactions vanish, 
we will end up with a new form of exactly solved Heisenberg model in which DM 
interactions are also present. The difficulty with this approach is that if we start with 
a uniform magnetic field, the final magnetic fields will not be uniform anymore, due 
to the periodic rotations. On the other hand if we start with a non-uniform magnetic 
field and demand that the final magnetic field will be uniform along the chain, then 
the spectrum of the original Hamiltonian will be very difficult to find. We hope to 
overcome these difficulties in our future work. 
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